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I.  Introduction 


Camouflage  paints  designed  to  aid  in  the  concealment  of 
military  objects  in  the  visible  and  near  infrared  spectral 
regions  would  only  fortuitously  accomplish  this  task  in  other 
regions  of  the  spectrum  as  well.  However,  it  has  recently 
become  a  major  concern  of  various  parts  of  the  United  States 
Department  of  Defense  and  specifically  the  MERADCOM  camouflage 
laboratory  to  be  able  to  develop  paints  or  coatings  that  would 
simultaneously  be  effective  in  the  thermal  infrared  region  of 
the  spectrum,  and  the  shorter  wave  near  infrared  portion  of 
the  spectrum  while  not  degrading  their  performance  in  the 
visible  region.  The  vital  importance  of  the  thermal  infrared 
region  which  we  will  define  here  as  principally  the  8-14  pm 
region  but  including  the  3-5  pm  region  as  well,  has  become  of 
continually  greater  concern  owing  to  the  rapid  development  of 
thermal  imaging  systems  notably  FLIRS,  and  other  line  scan¬ 
ners  in  which  the  platform  provides  the  extra  dimension  for 
construction  of  the  image.  The  extreme  sensitivity  of  these 
devices  and  their  present  availability  even  in  the  commercial 
world  has  dictated  that  scientists,  engineers,  and  paint  tech¬ 
nologists  concerned  with  these  problems  devote  a  considerable 
amount  of  effort  to  improving  the  camouflage  in  the  thermal 
infrared  region. 

This  problem  is  technologically  related  to  one  that  Arthur 
D.  Little  scientists  have  worked  on  over  the  last  decade  -  that 
of  infrared  remote  sensing  of  the  composition  of  a  planetary 
surface.  In  our  previous  work  on  the  problem  of  planetary 
surface  remote  compositional  sensing,  we  discovered  that  while 
suitable  data  in  the  infrared  region  could  readily  be  obtained, 
there  was  no  comprehensive  theoretical  method  available  to 
predict,  model  or  fully  understand  the  ways  in  which  the  signa¬ 
ture  of  a  mixture  of  materials  would  vary  with  the  amounts. 
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particle  sizes,  and  components  of  the  mixture.  To  fill  this 

need,  we  constructed  a  theoretical  model  of  reflectance  or 

emittance  of  a  mixture  of  minerals^  and  carried  out  a  number 

of  experimental  measurements  on  powders  to  aid  in  construe- 

2  3 

tion  of  the  model  '  .  The  only  existing  exact  theory  related 

4 

to  such  problems  is  the  Mie  theory  but  it  in  fact  only 
describes  the  absorbing  and  scattering  behavior  of  isolated 
spherical  particles.  Problems  rendering  this  theory  incor¬ 
rect  for  mixtures  of  mineral  particles  and  indeed  for  paints 
are  related  both  to  the  Mie  theory  requirements  for  spheri¬ 
city  and  for  isolation  of  the  individual  particles.  The 
theory  that  we  developed  treats  the  individual  particles  as 
being  randomly  oriented,  jagged  in  nature,  and  uses  a  modi¬ 
fied  six  beam  theory  of  radiative  transfer  related  to  the 

C 

Kubelka-Munk  formulation.  The  entire  theory  is  described 
in  a  number  of  publications  '  '  but  for  the  purposes  of  this 
report  it  is  sufficient  to  note  that  computer  calculations 
based  on  this  theory  successfully  simulate  to  a  high  degree 
the  spectral  signatures  of  various  minerals  and  mixtures  of 
minerals  including  lunar  soils®  such  that  we  can  now  success¬ 
fully  model  the  expected  spectrum  from  mixtures  of  materials 
of  varying  volume  fractions  and  particle  sizes  or  particle 
size  distributions. 

The  theory  is  made  up  of  two  subtheories,  one  for  parti¬ 
cles  large  compared  to  the  radiation  wavelength,  and  one  for 
particles  small  compared  to  the  radiation  wavelength.  The 
coarse  particle  theory  is  based  on  geometrical  optics  for 
spheres  but  modified  by  the  wave  optical  effects  of  edges  and 
other  asperities  on  the  surface  of  the  spheres  and  by  the 
proximity  of  neighboring  particles.  The  fine  particle  theory 
is  based  on  the  Lorentz-Lorenz  theory  of  dielectrics  and  treats 
particles  small  compared  to  the  radiation  wavelength  by  recog¬ 
nizing  that  the  radiation  necessarily  samples  "averaged”  opti¬ 
cal  properties  using  the  Lorentz-Lorenz  theory  to  calculate 
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the  average.  The  two  theories  are  merged  by  an  empirical 
bridging  relationship  which  has  been  found  to  provide  better 
agreement  with  experimental  results  in  the  intermediate  region 
than  either  theory  alone. 

The  successful  comparisons  of  simulations  and  experi- 

3 

mental  results  show  that  in  great  detail  the  physical  mech¬ 
anisms  of  scattering  and  absorption  by  particulate  media  are 
well  understood  and  represented  by  this  theory.  This  under¬ 
standing  led  us  to  believe  that  a  similar  theory  could  be 
constructed  for  fibrous  materials.  The  modified  theory  proved 
to  be  quite  able  to  predict  the  spectral  reflectance  or  emit- 
tance  of  fibrous  materials  as  a  function  of  fiber  diameter. 

7 

filling  factors  and  the  packing  in  individual  fiber  bundles  . 

These  two  successful  uses  of  the  general  theory  led  us  to 
believe  that  the  modeling  of  a  paint  or  coating  would  also 
succeed.  We  therefore  began  to  carry  out  the  necessary  modi¬ 
fications  of  the  theory  for  work  in  the  near  infrared  region 
of  the  spectrum  under  a  previous  contract.  At  MERADCOM's 
request,  we  enlarged  the  scope  of  this  previous  work  to  include 
the  thermal  infrared  region  of  the  spectrum  owing  to  the  vital 
nature  of  this  region  in  present  day  military  considerations. 

It  was  decided  that  we  should  establish  the  accuracy  of  the 
theory  and  refine  it  where  necessary  by  modeling  one  of  the 
Army's  existing  paint  formulations;  the  one  chosen  is  Forest 
Green.  However,  early  in  this  work  some  obvious  practical 
considerations  became  apparent  with  regard  to  the  nature  of 
the  problem  of  designing  improved  paints  or  coatings  for  this 
region.  These  will  be  discussed  later  in  the  body  of  the 
report. 

In  Section  II,  we  will  discuss  the  theoretical  changes 
required  to  convert  the  pre-existing  model  of  the  reflectance 
or  emittance  of  particulate  surfaces  to  the  case  of  particles 
in  a  binder  as  in  camouflage  paint.  In  Section  III,  we  will 
discuss  the  experimental  part  of  our  program  -  the  Forest  Green 
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paint  formulation  and  the  subformulations  that  we  obtained 
in  order  to  carry  out  some  modeling  prior  to  obtaining  all 
the  required  input  parameters  for  all  of  the  components  of 
Forest  Green  paint.  In  Section  IV,  we  present  the  optical 
constant  methodolbgy  and  results  during  our  first  six  months 
effort.  Finally,  in  Section  V,  we  will  show  the  fit  between 
theory  and  experiment  and  in  Section  VI,  present  suggestions 
for  further  work. 
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We  began  with  the  theory  as  presented  in  references  1,  2, 
and  3.  The  most  fundamental  theoretical  changes  required 
were  to  provide  for  the  occurence  of  surface  reflectance  terms 
originating  in  the  external  reflection  of  the  incident  beam 
by  the  Fresnel  relations.  It  has  always  been  part  of  our 
methodology  to  calculate  reflectance  rather  than  emittance 
as  it  is  intuitively  more  obvious  to  most  readers  as  well  as 
to  the  authors.  However,  the  emittance  can  be  obtained  from 
the  reflectance  by  the  use  of  Kirchhoff's  law  for  opaque 
substances.  That  law  which  is  that  absorptance  equals  emittance 
leads  to  the  formula? 

e  =  1-R 

Where  e  is  the  emittance  of  the  material  and  R  is  its  reflect¬ 
ance.  Both  quantities  may  be  subscripted  with  the  angular 
specifications.  That  is,  the  Kirchhoff  complement  of  the 
reflectance  at  some  angle  must  be  the  emittance  at  that  same 
angle.  As  reflectance  requires  two  subscripts  and  emittance 

g 

only  one,  analysis  shows  that  emittance  may  be  treated  as  if 
the  incident  beam  was  coming  from  inside  the  medium  and  entirely 
diffuse.  That  is,  the  proper  complement  of  the  emittance  at 
30°  is  the  reflectance  in  which  the  surface  is  illuminated 


diffusely  (hemispherically)  and  reflectance  is  measured  at 
30°  or  conversely  the  surface  is  illuminated  at  30°  and  the 
reflectance  is  measured  diffusely. 

The  surface  term  Re  in  our  present  formulation  refers  to 
the  value  obtained  by  integrating  the  Fresnel  relations  over 
all  angles,  i.e.,  diffuse  reflectance.  Associated  with  this 
diffuse  term  for  external  Fresnel  reflectance  is  an  internal 
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In  addition  to  these  surface  terms,  we  also  had  to  derive 
a  formula  including  the  volume  reflectance  resulting  from  a 
finite  paint  layer  on  a  substrate  of  reflectance  R  .  This 

g 

formula  is: 

R  -  IRe(l-RiRv)  (l-RvRg)-Re(VRi)  <W e'2TX 

+  d-Re>  (l-Ri)Rv(l-RvRg)-(l-Re)  (l-R^  (^,-Rg)  e'2yXl 
/  [  (1-R^RV)  U-RvR9)-(Rv-R1)  (Rv-Rg)e'2YX) 

where  x  is  the  thickness  of  the  paint  layer,  Ry  is  the  volume 

diffuse  reflectance  for  an  infinitely  thick  layer  and 
2  1/2 

y  =  (K  +  2KS)  '  *  K  and  S  are  the  usual  absorption  and  back- 

scattering  coefficients  given  in  our  previous  work1.  The 

derivation  of  this  formula  is  given  in  Appendix  1. 

Strictly  speaking,  R^  refers  to  the  reflectance  of  the 

underlying  substrate  as  modified  by  the  overcoating  layer. 

However,  as  the  reflectance  of  metal  surfaces  are  not  greatly 

affected  by  this  consideration,  we  have  used  the  directly 

observed  reflectance  (1-e)  of  our  substrate  aluminum  plate 

xn  our  calculations  thus  far.  The  R_  and  R.  for  the  fine 
,  e  i  i 

particle  theory  are  calculated  from  the  Lorentz-Lorenz 
complex  index  of  refraction  but  in  the  coarse  particle  theory 
we  have  made  the  assumption  that  the  only  coherent  surface 
terms  arise  from  the  matrix  so  that  these  terms  are  calculated 
for  the  matrix  material  as  if  it  were  totally  covering  the 
surfaces  of  the  pigment  particles.  All  pigment  contributions 
are  included  in  the  volume  term.  Finally,  as  our  previous 
theories  always  assumed  particle-vacuum  or  particle-air 
interfaces  only,  the  change  was  made  to  particle-matrix  inter¬ 
faces  by  dividing  particle  complex  indices  of  refraction  by 
the  real  part  of  the  matrix  complex  index  of  refraction. 
Wavelength  dependent  effects  were  likewise  adjusted  to  the 
index  in  the  medium  where  necessary. 
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III.  Experimental  Measurements 

As  previously  mentioned,  we  have  chosen  to  carry  out  our 
modeling  experiments  using  the  Army's  Forest  Green  paint  at 
the  suggestion  of  Dr.  Jay  Fox  then  of  MERADCOM.  This  is  a 
complex  material  composed  of  an  alkyd  binder  -  four  pigments 
and  two  extender  pigments.  The  composition  is  shown  in  Table  1. 


Table  I 


Compositions  of  Forest  Green  Formulations 


Component 


Volume  Fractions  in  Dry  Paint 


Forest  Green  (Sub  Formulations) 
MIL-E52798  FG2  FG3  FG4  FG6 

1)  Acid  Insoluble  Green 
Pigment  (Ferro  V-11655) 


(Co,Cr,Zn,Ti  Spinel) 

0.101 

0.103 

0.106 

0.117 

0.133 

2) 

Chrome  Yellow  (med) 

(PbCr04) 

0.026 

0.027 

- 

- 

- 

3) 

Molybdate  Orange 

0.008 

- 

- 

- 

- 

4) 

Cyandur  Violet  (Carbazole 
Dioxazine  Violet) 

0.010 

- 

- 

- 

- 

5) 

Talc  (Nytal  300) 

0.093 

0.093 

0.097 

- 

- 

6) 

Diatomaceous  Earth 
(Celite  110) 

0.175 

0.175 

0.183 

0.203 

- 

7) 

Alkyd  Resin  (Solids  Only) 
(Beckosol  11-070) 

0.587 

0.602 

0.614 

0.680 

0.867 
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Previous  work  under  joint  sponsorship  of  MERADCOM, 

NARADCOM  and  the  Army  Research  Office  indicated  that  the  near 
infrared  spectrum  of  this  material  is  dominated  by  the  acid 
insoluble  green  pigment,  as  is  true  to  a  lesser  extent  in  the 
visible  region.  In  the  visible  region  for  which  this  material 
was  originally  designed,  the  various  ingredients  performed 
their  functions  to  suitably  blend  with  the  chlorophyll  spectrum. 
If  it  were  only  for  the  near  infrared  region,  none  of  the  other 
tinting  pigments  would  be  required  and  only  the  extender 
pigments  whose  purpose  it  is  to  make  the  paint  into  the  matte 
finish  i.e.,  scatter  diffusely  would  be  required.  We  will 
digress  at  this  point  to  observe  that  the  WARSAW  pact  nations 
allow  their  paints  to  have  some  gloss  in  contrast  to  the  way 
our  paints  are  designed.  The  acid  insoluble  green  pigment 
(AIG)  is  a  mixed  oxide  of  zinc,  chromium,  titanium,  and  cobalt 
having  a  spinel  structure.  The  transmission  spectrum  of 
powdered  AIG  dispersed  in  a  pellet  made  of  thallium  chloride 
is  shown  in  Figure  1,  while  a  reflection  spectrum  obtained  from 
a  sample  which  was  hot  pressed  to  high  density  and  polished  is 
shown  in  Figure  2.  A  pellet  spectrum  of  lead  chromate  is  shown 
in  Figure  3  and  the  reflection  spectrum  of  a  single  crystal  of 
this  compound,  the  relatively  rare  mineral  known  as  crocoite, 
obtained  from  the  Harvard  University  collection  is  shown  in 
Figure  4.  Pellet  spectra  of  the  diatomaceous  earth,  the  talc, 
and  the  carbazole  dioxazine  violet  are  shown  in  Figures  5,  6  and 
7,  while  the  transmission  spectrum  of  a  thin  film  of  the  alkyd 
resin  is  shown  in  Figure  8.  A  spectrum  of  molybdate  orange 

9 

from  the  literature  is  shown  in  Figure  9. 

In  order  to  obtain  the  optical  constants  (complex  index 
of  refraction)  of  these  various  components  of  the  Forest  Green 
paint,  quantitative  data  had  to  be  obtained  and  reduced.  Our 
method  of  reduction  from  reflectance  spectra  is  given  in  refer¬ 
ence  10.  None  of  the  pellet  spectra  are  suitable  for  such  a 
purpose  except  with  large  relative  error.  Inasmuch  as  a  quan¬ 
titative  treatment  of  transmission  through  a  pellet  is  very 

8 


Arthur  D  Little  Inc 


A./um 


8 


H 

CC 

< 

U 

J) 

o 

o 

u 


•  9U»||IUI«U«UX 


13 


Arthur  DUttte  Inc 


4000  3600  3200  2800  2400  2000  *00  1200  800  490 

V,  CTO'* 

ALKYD  RESIN  FILM  TRANSMITTANCE  AT  BREWSTER  ANGLE 

Figure  8 
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difficult  owing  to  the  scattering  that  occurs  at  the  inter¬ 
faces  between  the  particles  and  the  medium  in  which  these  are 
dispersed.  As  will  be  seen  in  the  next  section,  the  reflec¬ 
tion  spectra  of  the  acid  insoluble  green  pigment  and  the 
crocoite  were  used  to  obtain  the  optical  constants  of  these 
materials  and  the  transmission  spectrum  of  the  alkyd  resin 
was  used  invoking  a  new  technique  (vide  infra)  to  obtain  the 
optical  constants  of  this  resin  binder. 

In  order  to  carry  out  our  modeling  calculations  in  the 
absence  of  detailed  refractive  indices  of  the  other  components 
and  also  to  gain  further  insight  into  the  properties  of  the 
Forest  Green  formulations,  we  asked  for  and  were  supplied  by 
MERADCOM  with  a  number  of  subformulations  (formulations  with 
various  ingredients  left  out'  of  the  Forest  Green  paint. 

These  include  what  we  refer  to  as  FG2,  FG3,  FG4 ,  and  FG6 ,  all 
listed  in  Table  I.  The  last  material  was  used  in  this  work 
for  a  first  simulation  in  the  Forest  Green  system.  It  should 
be  noted  that  as  FG6  has  no  extender  pigment,  it  is  rather 
glossy  appearing  in  the  visible  region.  However,  it  must  be 
noted  that  in  the  thermal  infrared  the  distinction  between 
pigments  and  extender  pigments  (those  which  do  not  contribute 
color)  is  completely  lost  as  both  the  talc  and  the  diatomaceous 
earth  have  strong  absorption  bands  and  so  produce  infrared 
"color." 

The  thermal  emittance  spectra  of  each  of  these  paints 

except  FG3  was  obtained  by  painting  a  layer  on  an  aluminum 

plate,  and  the  emittance  spectra  obtained  with  our  Michelson 

2 

Interferometer  Spectrometer  are  shown  in  Figure  10.  Another 
set  of  runs  on  FG4  at  a  greater  layer  thickness  (74  pm  vs  44  pm) 
showed  washed  out  spectral  features  with  a  slightly  higher 
emittance  level.  We  believe  the  washed  out  features  result 
simply  from  the  closeness  of  approach  to  an  emittance  of  unity 
caused  by  the  greater  layer  thickness.  Included  in  our  experi¬ 
mental  work  on  the  emittance  of  these  materials  in  the  7-35  pm 
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Forest  Green  (55/iin)  on  Aluminum 


10  900  TOO  too  1100  1300 

v,  cm'1 


EMITTANCE  SPECTRA  OF  FOREST  GREEN  SYSTEMS 

Figure  10 
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range  were  measurements  made  on  the  underlying  substrate 
aluminum  plate  and  measurements  made  on  the  alkyd  binder  alone 
and  a  possible  substitute  polyurethane  binder  (Figure  11) . 

All  of  these  spectra  were  obtained  in  duplicate  runs  with  the 
average  deviation  between  the  runs  being  0.05  (+0.025).  These 
spectra  were  obtained  using  two  blackbodies  at  different  temp¬ 
eratures  as  calibration  sources. 

The  samples  of  Forest  Green  paint  and  subformulations 
FG4  and  FG6  were  then  run  on  our  Digilab  Fourier  transform 
spectrometer  in  reflectance  and  the  spectra  are  shown  in 
Figures  12  through  14.  These  spectra  are  all  necessarily  qual¬ 
itative  inasmuch  as  no  reliably  known  reflectance  standard  was 
at  hand.  Some  of  these  were  run  against  the  reflectance  of 
the  unpolished  substrate  material  (aluminum)  and  some  against 
a  coarse  aluminum  scatter  plate  in  our  possession.  The  latter 
is  probably  a  somewhat  better  standard  as  it  is  a  more  diffuse 
reflector  having  an  almost  spectrally  flat  value  of  0.2  from 
about  600  to  4000  cm  1.  Clearly  a  diffuse  reflector  is  a 
requirement  for  a  standard  to  be  used  with  diffuse  samples  as 
a  specular  standard  does  not  allow  for  the  energy  reflected 
outside  of  the  acceptance  cone  of  the  optical  system.  The 
problem  then  becomes  one  of  relating  the  angular  distribution 
of  the  sample  and  standard  properly. 

As  sufficient  work  to  precisely  calibrate  our  reference 
plates  has  not  been  carried  out  we  consider  all  spectra  meas¬ 
ured  in  reflection  in  the  Digilab  instrument  to  be  qualitative 
in  nature  and  the  easiest  method  of  making  these  quantitative 
would  be  to  refer  them  by  Kirchhoff's  law  to  our  emittance 
spectra  in  the  7-35  ym  range.  The  greater  spectral  detail 
observed  in  reflectance  is  largely  the  result  of  scale  expan¬ 
sion  and  a  better  signal-to-noise  ratio. 

Comparisons  of  various  spectra  indicate  that  the  qualita¬ 
tive  appearance  of  the  various  bands  follows  Kirchhoff's  law 
well  (e.g..  Figure  15)/  and  that  the  data  are  quite  accurate 
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REFLECTANCE  OF  FOREST  GREEN  6  VS  DIFFUSE  REFERENCE 


Arthur  D  Little  Inc 


COMPARISON  OF  EMITTANCE  AND  REFLECTANCE  DATA 

FOR  FOREST  GREEN  6 
Figure  15 
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with  respect  to  spectral  shape  over  the  whole  range  of  the 
measured  spectrum. 

Reflectance  data  from  the  alkyd  binder  is  shown  in 
Figure  16.  The  low  reflectance  shown  by  this  film  in  the 
8-14  urn  region  should  be  contrasted  with  the  relatively 
higher  value  in  the  3.5-5. 5  urn  region.  In  the  same  figure, 
we  show  a  comparable  spectrum  of  the  polyurethane  binder 
that  has  been  suggested  as  a  possible  replacement  for  the 
alkyd  binder  for  contamination  resistance11.  The  data  on 
the  polyurethane  binder  was  obtained  after  it  became  apparent 
from  our  spectra  that  no  combination  of  pigments  could  sub¬ 
stantially  reduce  the  thermal  emittance  in  the  8-14  urn  region 
below  approximately  0.9  inasmuch  as  the  alkyd  binder  itself 
had  that  level  of  emittance  in  this  region.  As  can  be  seen 
from  Figure  16,  the  polyurethane  binder  is  no  better  and  as 
a  consequence  of  these  results  we  are  led  to  state  that  no 
paint  will  provide  a  suitably  low  emittance  in  the  8-14  pm 
region  if  it  is  produced  using  either  the  alkyd  or  polyure¬ 
thane  binders. 

In  fact,  it  may  be  very  difficult  to  find  an  appropriate 
organic  binder  if,  as  we  believe,  it  is  necessary  to  reduce 
the  emittance  from  values  close  to  unity  to  values  of  less 
than  0.3  over  the  appropriate  thermal  band.  We  have  suggested 
that  measurements  be  made  in  the  future  on  paraffins  to 
ascertain  whether  these  are  feasible  substitute  binders,  in 
the  sense  of  being  capable  of  providing  a  suitable  physical 
material  for  a  paint  binder  and,  if  that  is  possible,  of 
providing  a  sufficiently  low  binder  absorption  such  that  the 
binder  itself  will  not  constitute  a  primary  source  of  thermal 
emission.  A  simple  way  to  view  this  problem  is  that  by 
Kirchhoff's  law,  the  emittance  is  equal  to  the  absorptance 
so  that  any  substance  that  absorbs  well  will  emit  well.  Thus 
the  problem  is  reduced  to  one  of  finding  a  sufficiently  trans¬ 
parent  binder  that  the  binder  itself  does  not  constitute  a 
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Figure  16 


thermally  bright  source  and  thus  may  be  used  as  a  matrix  to 
contain  materials  which  scatter  the  radiation  sufficiently 
to  allow  for  a  low  emittance  composite  material.  We  will 
discuss  this  again  further  in  Section  VI. 
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IV.  Optical  Constants  (Complex  Refractive  Indices) 


In  order  to  carry  out  our  theoretical  simulations ,  we 
require  as  input  parameters 

1.  particle  size  or  size  distribution 

?.  volume  fractions  of  the  various  components,  and 

3.  the  spectral  complex  refractive  index  of  each 
material. 

We  were  able  to  obtain  such  information  during  this  contract 
by  using  our  previously  reported  technique10  for  a  hot 
pressed  sample  of  the  acid  insoluble  green  pigment  (fit 
shown  in  Figure  2) ,  and  the  single  crystal  crocoite  (fit 
shown  in  Figure  4)  sample.  As  the  molybdate  orange  and 
carbazole  dioxazine  violet  are  present  only  in  very  small 
quantities,  we  chose,  owing  to  the  present  limitations  of 
time  and  funding  to  neglect  these  for  our  initial  compar¬ 
isons.  We  attempted  to  obtain  a  talc  sample  with  sufficient 
polish  to  obtain  a  good  reflection  spectrum  but  some  diffi¬ 
culty  was  obtained  in  getting  a  sufficiently  specular  surface 
and  so  that  material  had  to  be  abandoned  for  the  time  being. 
Qualitative  measurements  based  on  the  pellet  spectra  could 
of  course  be  used  to  get  approximate  optical  constants  for 
any  of  these  materials,  but  we  chose  to  try  to  obtain  the 
best  optical  constants  possible  for  the  first  paint  simula¬ 
tions.  The  diatomaceous  earth  presents  an  extreme  problem 
in  obtaining  a  homogeneous  smooth  sample  for  optical  constant 
determination.  For  the  relatively  transparent  alkyd  binder, 
a  reflection  spectrum  does  not  show  sufficiently  strong  bands 
to  be  able  to  use  our  previous  technique  of  fitting  specular 
reflectance  to  obtain  optical  constants.  Therefore,  we  chose 
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to  use  a  transmission  spectrum  of  a  thin  film  of  this  mater¬ 
ial.  Some  difficulty  was  encountered  in  attempting  to  obtain 
a  sufficiently  smooth,  homogeneous  film  of  the  binder.  Reduc¬ 
tion  of  the  data  from  our  films  by  the  usual  transmission 
techniques  involving  the  full  transmittance  equation,  includ¬ 
ing  reflection  losses,  was  considered  likely  to  involve  some 
error,  owing  to  the  lack  of  optical  perfection  of  the  film 
surfaces. 

In  order  to  circumvent  this  difficulty,  we  have  chosen  to 
use  a  somewhat  different  technique  based  on  measurements  of 
the  transmission  of  the  best  quality  film  that  we  could  obtain, 
at  the  Brewster  angle.  The  Brewster  angle,  <j>  is  defined  by 
n  *  tan  $,  where  n  is  the  real  part  of  the  complex  refractive 
index  of  the  medium.  At  the  angle  of  incidence  $ ,  radiation 
whose  electric  vector  is  in  the  plane  of  incidence  cannot  be 
reflected,  so  for  this  particular  angle,  and  with  the  use  of 
a  polarizer  that  passes  only  such  radiation,  the  reflection 
loss  will  be  zero.  That  being  true,  the  usual  law  for  thick 
samples  (without  interference  phenomena) 

T  =  (1-R)  V°X 

l-R2e-2ax 

is  simplified  to  T  *  e~aX.  Thus  we  were  able  to  avoid 
problems  involving  a  least  squares  fit  to  the  reflection 
losses  where  the  reflectance  is  from  an  imperfect  surface. 

Of  course  no  polarizer  is  perfect  and  spectral  variations  in 
the  value  of  n  would  in  principle  lead  to  varying  values  of 
Brewster' 8  angle.  For  organic  films,  however,  the  disper¬ 
sion  of  the  refractive  index  is  quite  small  and  an  error 
analysis  indicated  that  the  method  should  be  quite  accurate. 

The  results  of  our  standard  technique  for  fitting  reflec¬ 
tance  data  to  a  set  of  Lorentz  line  parameters  (classical 
dispersion  theory  oscillators)  are  shown  in  Figures  2  and  4 
for  the  acid  insoluble  green  hot  pressed  sample  and  the 
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crocoite  sample  measured  at  30°  angle  of  incidence.  The  crystal 
orientation  for  the  latter  was  9°  from  the  100  face.  This  was 
the  only  face  we  could  obtain  for  this  monoclinic  crystal  in 
sufficiently  large  size  for  our  measurements.  An  attempt  to 
hot  press  a  sample  of  lead  chromate  powder  in  the  same  way 
as  carried  out  with  the  acid  insoluble  green  pigment  failed 
owing  to  chemical  reactivity.  Therefore,  the  Lorentz  line 
parameters  found  from  this  arbitrary  face  of  the  crocoite 
sample  represent  some  combination  of  the  three  sets  of  line 
parameters  that  should  be  used  to  represent  this  monoclinic 
material.  However,  the  principal  line,  which  is  found  to  be 
at  823  cm”*,  Table  II,  is  seen  to  correspond  roughly  to  the 
observed  transmission  of  the  pellet  spectrum  (Figure  3)  which 
shows  one  large  absorption  band  centered  near  840  cm  *.  The 
figure,  however,  does  indicate  some  evidence  of  structure  in 
this  band  which  probably  indicates  contributions  from  the 
other  orientations  of  the  material.  A  lower  intensity  feature 
near  400  cm”*  did  not  provide  a  sufficiently  strong  band  in 
the  reflection  spectrum  of  our  orientation  of  crocoite  to 
utilize  a  second  line  in  our  Lorentz  line  parameter  fit 
although  there  is  some  evidence  for  its  presence  (see  Figure  4) . 
A  fit  with  more  lines  could  probably  be  attempted  but  was  not 
carried  out  as  crocoite  was  not  used  in  our  initial  modeling. 

The  Lorentz  line  fit  to  the  acid  insoluble  green  sample 
is  shown  in  Figure  2,  and  the  Lorentz  line  parameters  are 
given  in  Table  III.  The  strength  of  the  two  lines  will  clearly 
dominate  the  model  spectrum  in  the  low  frequency  region  so 
that  the  small  crocoite  line  mentioned  above  is  undoubtedly 
unnecessary. 

For  the  alkyd  resin,  a  thin  film  having  a  thickness  of 
2.6  ym  was  fitted  by  the  Brewster  angle  technique.  The  Lorentz 
line  parameters  for  the  alkyd  film  are  given  in  Table  IV,  and 
the  fit  is  shown  in  Figure  17.  As  can  readily  be  seen,  the 
strengths  of  these  lines  are  much  less  than  those  of  the  acid 
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insoluble  green  and  the  crocoite  samples.  This  is  the  reason 
that  insufficiently  strong  reflectance  spectra  could  be 
obtained,  but  the  transmission  technique  can  be  seen  to  have 
worked  quite  well.  As  will  be  shown  later,  in  the  thicknesses 
of  paint  layers  commonly  used,  these  lines  are  sufficiently 
strong  to  produce  features  and  to  emit  strongly. 


Table  II 

Crocoite 

Estimated  standard  deviations  are  given  in  parentheses 

e  =  6.2066  (0.0646) 

00 

v  S  Y 

823.32  (1.53)  1.68338  (0.03449)  0.05146  (0.00183) 


Table  III 


Acid  Insoluble  Green  Pigment 


Estimated  standard  deviations  are  given  in  parentheses 

e  =  4.5228  (0.0800) 

00 


491.93  (1.04) 

611.37  (0.79) 


3.60492  (0.06646)  0.03733  (0.00219) 
0.66597  (0.02100)  0.03409  (0.00163) 
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The  Brewster  angle  technique  which  we  have  used  in  this 
is,  we  believe,  novel  and  will  adequately  produce 
Lorentz  line  parameter  fits  for  relatively  weakly  absorbing 
materials  in  the  same  manner  that  the  stronger  materials  can 
be  fitted  by  reflection  techniques.  We,  therefore,  believe 
that  we  have  techniques  for  almost  any  kind  of  sample  if  a 
homogeneous  film  or  plate  can  be  prepared. 
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V.  Comparison  of  Theoretical 
Simulations  with  Experimental  Measurements 


As  a  test  case  of  the  ability  of  our  model  to  predict 
the  spectral  signature  of  the  Forest  Green  paint  system,  we 
have  chosen  the  FG6  paint,  which  consists  only  of  the  alkyd 
binder  and  the  acid  insoluble  green  pigment  for  modeling 
purposes.  Figure  18  shows  the  particle  size  distribution  of 
the  acid  insoluble  green  pigment  supplied  to  us  by  Mr.  Aladar 
Burgyan  of  the  Ferro  Corporation.  With  this  curve  as  a  guide, 
we  selected  monodisperse  particle  diameters  of  1,  2,  and  5  urn 
for  the  purpose  of  modeling  the  variations  with  particle  size 
that  may  be  expected  in  the  range  of  the  particle  sizes  pre¬ 
sent.  The  results  of  our  first  simulations  are  shown  in 
Figure  19  for  these  three  particle  sizes.  These  may  be  com¬ 
pared  with  the  emittance  measurements  on  the  FG6  formulation 
made  by  our  Michelson  interferometer  and  the  qualitative 
reflectance  measurements  made  on  our  Digilab  system  shown  in 
Figure  15. 

The  emittance  data  are  redrawn  in  Figure  19  to  facilitate 
the  detailed  comparison.  The  important  point  to  note  is  the 
approximate  match  in  overall  spectral  level  and  the  generally 
accurate  major  features  near  500  and  630  cm-1  resulting  from 
the  AIG  pigment.  The  discrepancy  in  spectral  level  observed 
(AR=Ae''-0.1  or  10%  in  emittance)  is  partially  due  to  the  use  of 
R  =  0.9  (From  Figure  11)  rather  than  the  value  for  the  over- 
coated  substrate.  It  also  is  likely  to  involve  some  of  our 
experimental  error,  as  all  the  emittance  spectra  shown  in  this 
report  were  from  typical  single  runs  rather  than  averaged  with 
the  duplicate  runs.  That  some  of  the  detailed  small  discrep¬ 
ancies  observed  are  probably  the  result  of  small  inaccuracies 
in  our  bridging  relationships  is  shown  in  Figure  20  where  the 
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separate  contributions  of  the  coarse  and  fine  particle 
theories  are  shown  together  with  the  combined  result.  It 
can  be  seen  that  the  fine  particle  theory  produces  the  best 
qualitative  fit  to  the  location  and  relative  importance  of 
the  various  spectral  features  but  the  poorest  level  when 
compared  to  the  emittance  data.  A  somewhat  modified 
bridging  rule  would  be  expected  to  improve  the  overall 
situation. 
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VI 


Suggestions  for  Further  Work 


As  discussed  above ,  we  consider  the  basic  problem  of 
reducing  the  thermal  emittance  of  any  paint  to  be  principally 
one  of  finding  a  relatively  transparent  binder  in  the  thermal 
infrared  region.  As  existing  binders  tend  to  be  quite  opaque 
even  in  very  thin  layers,  we  have  discussed  with  Mr.  Buccigross 
of  our  staff,  the  possibility  of  finding  a  new  binder  that 
might  satisfy  the  transmission  requirement.  He  posed  for  us 
a  number  of  possible  materials,  and  after  having  surveyed  his 
list  and  also  the  infrared  characteristics  expected  for  such 
films,  we  have  come  to  the  conclusion  that  the  best  organic 
type  binder  would  necessarily  be  one  involving  some  type  of 
paraffin,  polyethylene,  or  rubber  material.  These  materials 
are  not  presently  available  in  a  form  suitable  fo^  making  a 
paint,  and  so  in  order  to  develop  such  a  paint  or  a  coating 
some  development  work  along  the  lines  of  producing  such  a 
vehicle  would  be  required.  Inorganic  binders  should  also  be 
considered. 

Assuming  an  effective  bland  binder  could  be  developed,  it 
is  expected  that  the  visible  and  near  infrared  properties  would 
be  quite  similar  to  those  presently  obtained  with  the  existing 
camouflage  paints  by  using  the  same  pigments  especially  if  the 
binder  refractive  index  is  near  1.5.  In  the  thermal  infrared, 
however,  the  pigments  would  have  to  be  modified  as  the  present 
ones  are  far  from  ideal.  This  requirement  would  have  to  be 
satisfied  without  substantially  degrading  the  performance  of 
the  paints  in  the  visible  and  near  infrared  regions. 

Further  modeling  work  along  the  lines  begun  in  this  report 
should  be  carried  out  using  the  candidate  binders  and  improved 
pigments  required  for  the  improved  paint  formulations  in  order 
to  optimize  particle  sizes  and  volume  fractions. 
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OF  A  LAYER  OF  PAINT  ON  A  SUBSTRATE 


For  mathematical  convenience  we  consider  that  the  layer 
of  paint  consists  of  one  or  more  kinds  of  pigment  particles 
embedded  in  a  continuous  binder  which  extends  for  a  small 
distance  in  front  of  the  pigment  particles  and  also  behind 
them.  The  thin  film  of  pure  binder  at  the  front  gives  a 
well-defined  gloss  component  of  the  reflectance.  The  thin 
film  of  binder  at  the  back  is  assumed  to  be  in  good  optical 
contact  with  the  substrate  and  therefore  provides  a  definite 
value  of  the  back  reflectance  Rg.  We  also  assume,  again  for 
mathematical  tractability ,  that  the  light  is  everywhere 
completely  diffuse. 

The  radiative  transfer  in  the  paint  layer  can  be  described 
in  terms  of  the  well-known  two-flux  model  of  Kubelka  and  Munk^5^. 
This  theory  yields  the  general  solution: 

I  =  Ae“YX  +  BeyX  (Al) 

J  =  RvAe"YX  +  ^"BeYX  (A2) 

where  I  and  J  are  respectively  the  ingoing  and  outgoing  diffuse 
fluxes  at  a  distance  x  from  the  front  surface  of  the  paint 
layer. 

A  and  B  are  arbitrary  constants  which  are  determined  by 
the  boundary  conditions  at  the  front  and  back  of  the  layer. 
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Ry  is  the  volume  reflectance  of  the  pigment-binder  medium, 
i.e.,  the  ratio  J/I  of  the  outgoing  to  ingoing  flux  that  would 
exist  if  the  paint  layer  extended  to  infinity. 

The  quantity  y  is  the  extinction  coefficient  of  radiation 
in  an  infinite  medium,  and  includes  the  effect  of  both  scatter¬ 
ing  and  absorption  losses  in  the  medium. 

Rv  and  y  are  related  to  the  diffuse  absorption  and  scatter¬ 
ing  coefficients  K  and  S  by  the  relations: 


1/2 


(A3) 


y  =  K2  +  2KS 


(A4) 


For  given  values  of  the  complex  index  of  refraction 
m  »  n-ik  of  the  pigment  particles  and  the  binder,  and  the 
average  particle  size  and  packing  density,  we  can  calculate 
K  and  S  by  means  of  a  theory  '  '  developed  by  us  for  powders. 
We  have  to  modify  this  theory  for  the  case  of  particles 
immersed  in  a  binder  by  replacing  the  complex  index  m  of  the 
particles  by  the  relative  index  m/nD  of  a  particle  with  respect 
to  the  binder,  and  at  the  same  time  use  the  radiation  wave¬ 
length  X/nQ  in  the  binder  in  place  of  the  free-space  wave¬ 
length  X. 

The  boundary  conditions  needed  to  determine  the  constants 
A  and  B  in  Equations  (Al)  and  (A2)  will  now  be  considered.  At 
the  back  interface,  located  at  x  «  X,  the  boundary  condition 
is  simply: 

J (X)  «  Rg  I(X)  (A5) 
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At  the  front  surface  the  boundary  conditions  are  some¬ 
what  more  complicated.  The  incident  diffuse  radiation  strikes 
the  binder  surface  and  is  partially  reflected.  We  calculate 
the  corresponding  external  diffuse  reflectance  Rg  by  integrat¬ 
ing  the  Fresnel  reflectance  over  all  angles  of  incidence. 

The  fraction  of  the  radiation  entering  the  binder  is  therefore 

i-V 

Back- scattered  diffuse  radiation  arriving  at  the  binder 

surface  internally  is  also  partially  reflected.  The  internal 

diffuse  reflectance  R^  is  not  equal  to  Rg.  If  a  piece  of 

dielectric  material  of  index  nQ  is  imagined  to  be  put  in  a 

cavity  at  some  temperature  T  the  radiation  flux  in  the  dielec- 

2 

trie,  under  thermodynamic  equilibrium  conditions,  is  nQ 
times  as  large  as  in  the  empty  part  of  the  cavity.  This 
implies  that,  for  equilibrium  to  be  maintained,  the  external 
and  internal  transmission  coefficients  T_  and  T.  through  the 
surface  of  the  dielectric  must  satisfy  the  condition: 


T 


e 


It  follows  that: 


Ri 


(A6) 


(A7) 


In  terms  of  Re  and  R^  the  boundary  conditions  at  the 
front  surface,  x  =  o,  are: 


I(o)  =  d-Rg)  IQ  +  R^lo) 


(A8) 


Jo  -  Re  Io  +  U-*i>J<o) 


(A9) 
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where  I  and  J  are  the  incident  and  reflected  diffuse  fluxes, 
o  o 

The  diffuse  reflectance  R  is  given  by: 


R  = 


(A10) 


On  applying  the  above  boundary  conditions  to  Equations 
(Al)  and  (A2)  we  find#  after  some  algebra#  that: 


-2yx 


Re(1-RiV<l-VV  '  WRi)(VVe 

[4-  (l-Re)g-R1)Rv(l-RvRa)  -  (1-Re)  (1-^)  (R,-Rq)e  *"21 


-2yx 


<l-*iV  U-W  -  <VRi>(W« 


-2yx 


(All) 
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